Document  Number 
TRA.COR  66 -136 -U 


CO 

<X) 


TECHNICAL  MEMOtANDUM 


PARAMETRIC  ANALYSIS  OF  DOME-TRANSDUCER  INTERACTIONS 


Prepared  for  the 


Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington,  D.  C. 

Attention:  Code  1622E 


Contract  NObsr-91229 

Project  Serial  Number  SFlOl-OB-l? 

Task  8139 


TRACOR,  INC. 


6500  TRACOR  LANE.  AUSTIN.  TEXAS  78721 


Document  Number 
TRACOR  66- 136 -U 


TECHNICAL  MEMCilANDllM 


A  PARAMETRIC  ANALYSIS  OF  DOME -TRANSDUCER  INTERACTIONS 


Prepared  for  the 


Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington,  D.  C. 

Attention:  Code  1622E 


Contract  NObsr-91229 

Project  Serial  Number  SFIOI-O3-I7 

Task  8139 


February  10,  1966 


Prepared  by: 

Moyef 
Morel  1 


Document  Number 
TRACOR  66-136-U 


TRACOR,  INC.  6500  TRACOR  LANE,  AUSTIN,  TEXAS  78721 

ABSTRACT 

This  technical  memorandum  presents  some  results  of  a 
parametric  analysis  of  the  interaction  of  sonar  domes  and  trans¬ 
ducers  during  transmission.  The  results  were  obtained  from  the 
numerical  evaluation  of  an  analytical  model  consisting  of  a 
cylindrical  transducer  and  concentric  dome.  The  dome -transducer 
parameters  considered  in  this  study  include  dome -transducer  spac¬ 
ing,  transducer  diameter,  transmitting  frequency,  dome  thickness, 
and  material  properties  of  the  dome.  The  parameter  ranges  include 
values  appropriate  to  the  AN/SQS-4,  AN/SQS-23,  and  AN/SQS-26  systems. 
The  results  are  presented  as  functions  of  the  dimensionless  products 
of  rhe  dome -transducer  parameters.  In  this  form,  the  results 
represent  an  initial  effort  at  formulating  guidelines  which  have 
broad  application  in  dome  design.  Improved  guidelines  will  be 
generated  as  better  analytical  models  of  the  dome-transducer 
are  developed. 

The  results  show  that  the  dome  produces  a  degradation  of 
each  of  the  performance  criteria  compared  to  a  bare  transducer. 

The  effects  of  the  dome  can  be  reduced,  however,  by  proper  choice 
of  dome  parameters.  The  results  provide  a  quantitative  measure 
of  the  relative  importance  of  the  various  parameters  in  determin¬ 
ing  degradation  in  performance  due  to  the  dome.  It  is  shown 
that  the  product  of  dome  material  density  and  thickness  (surface 
density)  is  the  most  important  parameter  in  determining  the  effect 
of  the  dome  on  transducer  performance. 
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I.  INTRODUCTION 

A  set  of  design  goals  or  performance  criteria  and  an 
analytical  model  of  the  dome -transducer  are  two  of  the  ingredi¬ 
ents  necessary  in  the  development  of  guidelines  for  dome  design. 
Guidelines  are  generated  by  determining  the  functional  dependence 
of  each  of  the  performance  criteria  on  the  parameters  of  the  dome- 
transducer  system.  By  choosing  those  parameter  values  which 
optimize  the  performance  criteria,  the  designer  is  led,  hopefully, 
to  an  optimum  dome  design.  The  current  inadequacy  of  design 
information  is  due  basically  to  inadequate  models  for  dome- 
transducer  interactions. 

This  technical  memorandum  presents  the  results  of  an 
initial  effort  to  generate  dome  design  guidelines.  While  the 
dome -transducer  model  is  not  yet  completely  realistic  in  terms 
of  dome  geometry,  the  type  of  information  and  method  of  presenta¬ 
tion  of  the  results  should  be  representative  of  design  guideline 
requirements.  The  guidelines  presented  herein  were  generated  by 
computing  the  dependence  of  certain  performance  criteria  on  the 
parameters  of  a  dome -transducer  model  consisting  of  a  cylindrical 
transducer  and  concentric  dome.  The  performance  criteria 
considered  in  this  study  include: 

1.  Directivity  Index 

2.  Major  Lobe  Width 

3.  Side  Lobe  Level 

4.  Dome  Transmission  Coefficient 

5.  Dome  Transfer  Coefficient 

6.  Nearfield  Roughness. 

The  first  three  quantities  are  familiar  measures  of  the  farfield 
characteristics  of  a  dome -transducer. 

Dome  Transmission  Coefficient  is  a  measure  of  the  per¬ 
centage  of  energy  incident  on  the  dome  which  is  transmitted  to 
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the  far field.  A  portion  of  the  incident  energy  is  reflected 
within  the  dome  and  is  not  tra.ismitted  to  the  farfield.  These 
reflections  form  standing  waves  between  the  dome  and  transducer 
and  affect  the  coupling  of  tha  transducer  to  the  medium. 

Dome  Transfer  Coefficienc  is  a  measure  of  the  coupling 
between  transducer  and  dome  and  is  a  function  of  system  geometry. 
Proper  coupling  results  in  an  efficient  power  transfer  from  the 
transducer  to  the  medium. 

Nearfield  Roughness  is  related  to  the  transducer  source 
level.  In  a  typical  transducer,  element  mutual  interactions 
and  dome -transducer  interactions  produce  "roughness"  or  non¬ 
uniformity  in  the  nearfield  sound  pressure  distribution.  It  is 
desirable  to  maintain  a  smooth  nearfield  pressure  distribution; 
otherwise,  it  may  be  necessary  to  reduce  the  transducer  source 
level  in  order  to  keep  the  near field  sound  pressure  peaks  below 
the  cavitation  limit.  Consequently,  nearfield  rouglmess  can  be 
an  important  performance  criterion  consistent  with  obtaining 
msaximum  source  level. 

The  computed  results  for  performance  criteria  are  presented 
as  functions  of  dimensionless  products  of  the  dome -transducer 
parameters.  This  method  of  presentation  does  not  limit  the  results 
to  a  specific  sonar  system  but  enables  the  results  to  be  applied 
to  any  system  whose  configuration  is  compatible  with  the  parameter 
ranges  used  in  the  computations.  The  parameter  ranges  include 
values  appropriate  to  the  AN/SQS-4,  AN/SQS-23,  and  AN/SQS-26  systems. 

In  addition  to  the  specific  results,  some  general  conclusions 
concerning  the  effects  of  domes  on  sonar  performance  are  included. 

A  summary  of  the  results  and  their  implications  in  dome  design 
are  given  also. 
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II.  ANALYSIS 

This  section  contains  a  brief  discussion  of  the  dome- 
transducer  model  used  for  the  analysis  described  in  this  memorandum, 
the  definitions  of  the  performance  criteria,  and  the  formulation 
of  a  method  of  presentation  of  the  results. 

A.  Dome -Transducer  Model 

The  dome -transducer  model  consists  of  a  cylindrical 
transducer  of  radius  "a"  surrounded  by  a  homogeneous,  elastic 
shell  of  radius  "b”  and  thickness  "h.'*  (See  Fig.  1.)  The  problem, 
as  formulated,  is  two-dimensional.  Twenty-four  staves  of  the 
transducer  are  active  (radiating) ;  the  remaining  forty-eight 
staves  are  rigid.  The  active  staves  are  phased  to  a  straight  line 
and  have  a  cosine  amplitude  shading.  The  stave  width  is  approxi¬ 
mately  one-half  the  wavelength  at  the  transmitting  frequency. 

The  equations  describing  the  sound  pressure  field  inside  and 
outside  the  dome  as  a  result  of  the  radiation  of  the  transducer 
are  given  in  Appendix  A.  The  complete  development  of  the  analyti¬ 
cal  model  can  be  found  in  Ref.  1.  The  equations  given  in  Appendix 
A  are  used  to  compute  the  various  performance  criteria  studied 
in  this  memorandum. 

® •  Performance  Criteria 

The  dome -transducer  performance  criteria  which  are 
considered  in  this  memorandum  are  defined  as  follows: 

1.  Directivity  Index 

A  two-dimensional  Directivity  Index  can  be  defined 
by  the  expression 


Directivity  Index  =  20  log-j^Q 


p 

*^max 

P  ’ 
‘^avg 


(1) 
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FIG.  I  -  GEOMETRY  FOR  A  CYLINDRICAL  TRANSDUCER 
RADIATING  IN  A  CONCENTRIC  DOME 
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where  p  is  the  maximum  sound  pressure  in  the  farfield,  and 
^max  ’ 

^avg  average  sound  pressure  taken  over  360°  azimuth. 

The  Directivity  Index  is  computed  directly  from  the  numerical 
results  for  the  farfield  sound  pressure  field  of  the  dome-transducer . 


2..  Major  Lobe  Width 

The  Major  Lobe  Width  is  defined  as  the  width  of 
the  major  lobe  at  the  6  dB  down  points  of  the  farfield  directivity 
pattern.  Major  Lobe  Width  is  obtained  directly  from  the  numerical 
results  for  the  farfield  directivity  pattern. 


3.  Side  Lobe  Level 


Side  Lobe  Level  is  defined  by  the  expression 


Side  Lobe  Level  =  20  login 

^max 

where  is  the  sound  pressure  amplitude  of  the  highest  side 
lobe  in  the  farfield  directivity  pattern,  and  is  the  maximum 

sound  pressure  amplitude  in  the  farfield  directivity  pattern. 

4.  Dome  Transmission  Coefficient 

Dome  Transmission  Coefficient  is  given  by  the 

expression 


t 

Dome  Transmission  Coefficient  =  ^  .  (3) 

^i 

where  is  the  time-average  power  transmitted  through  the  dome, 
and  is  the  time -average  power  incident  on  the  dome.  The  power 
incident  and  the  pov/er  transmitted  are  determined  by  the  relations 
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2tt 

P  =  Re  f  p  (r,0)|  •  v  *(r,0)|  bd0 

^  "o  ^  r=b  ^  r=b 

and 


Pi(r, 0) I 

^  r=b 


v.*(r,0)j  bd0 
^  r=b 


(4) 


(5) 


where  p^(r,0)  and  v^*(r,0)  are  respectively  the  sound  pressure 
and  the  complex  conjugate  of  the  normal  particle  velocity  of 
the  sound  field  transmitted  through  the  dome,  and  P£(r,0)  and 
v^*(r,0)  are  the  sound  pressure  and  the  complex  conjugate  of 
the  normal  particle  velocity  of  the  sound  field  incident  on 
the  dome.  These  quantities  (p. ,  v  p.  and  v.*)  are  computed 
as  described  in  Appendix  B.  Dome  Transmission  Coefficient,  as 
defined  here,  is  primarily  a  function  of  the  dome  thickness  and 
material  properties  and  not  the  dome -transducer  spacing.  It  does 
not  measure  any  loss  in  signal  level  due  to  dome -transducer 
interaction,  but  simply  the  ratio  of  power  transmitted  through 
the  dome  to  power  incident  on  the  dome. 

5.  Dome  Transfer  Coefficient 


Dome  Transfer  Coefficient  is  defined  by  the 


expression 


'p 

Dome  Transfer  Coefficient  =  20  log^^Q  (  dome 


where  P 


with  dome 


no  dome 

t  is  the  time -average  power  transferred  to  the 


(6) 


^^with  dome  equal  to  P^  as  defined  for  Dome  Transmission 

Coefficient.  P.  is  of  arbitrary  level.  ^^j[th  dome  determined 
from  a  knowledge  of  the  transducer  element  velocities  and  the 
radiation  loading  at  the  transducer  elements  including  reflections 
from  the  dome. 
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water  by  the  transducer  with  dome,  and  P  j  ^  is  the  time- 
overage  power  transferred  to  the  water  by  a  bare  transducer. 

If  it  is  assumed  that  element  velocity  control  ercists  in  the 
transducer,  i.e.,  the  head  velocities  of  the  transducer  elements 
are  independent  of  radiation  impedance,  it  is  permissible  to 
compare  the  power  transferred  to  the  water  with  and  without  the 
dome.  (If  velocity  control  does  not  exist,  the  analysis  becomes 
much  more  complex,  the  head  velocity  for  each  element  must  be 
obtained  by  solving  a  large  set  of  simultaneous  equations. 

Absence  of  velocity  control  can  effect  a  degradation  of  sonar  per¬ 
formance.)  The  power  P  transferred  is  given  by  the  expression 

TT 

P  -  Re  J  p(r,0)|^,^  •  v’^(r,0)l^,g  adfi  ,  (7) 

-TT 

whore  p(r,0)  is  the  sound  pressure,  and  v  (r,9)  is  the  complex 
conjugate  of  the  normal  particle  velocity.  Both  p(r,9)  and 
V  (r,9)  are  obtained  from  the  expressions  given  in  Appendix  A.- 
The  integral  is  taken  over  the  transducer  face.  (See  Fig.  1.) 

Dome  Transrai’ ssion  Coefficient  measures  that  fraction  of  Incident 
energy  which  is  transmitted  through  the  dome.  The  remainder  of 
the  energy  is  reflected  from  the  dome.  Dome  Transfer  Coefficient 
is  a  measure  of  the  manner  in  which  the  reflected  energy  Interferes 
with  the  transducer. 

6.  Nearfield  Roughness 

Nearfleld  Roughness  is  defined  by  the  expression 

p 

Nearfleld  Roughness  -  20  login  ,  (8) 

p 

*^mean 

where  is  the  maximum  pressure  amplitude,  and  is  the 

mean  pressure  amplitude  on  the  active  portion  of  the  transducer 
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face.  Nearfield  Roughness  is  obtained  directly  from  the  numerical 
evaluation  of  the  nearfield  sound  pressure  of  the  dome -transducer 
model . 

C .  Parameters  of  the  Dome -Transducer  Model 

The  parameters  of  the  dome -transducer  model  and  their 
dimensions  are  given  in  the  table  below. 


PARAMETER 
Transducer  Radius 
Dome  Radius 
Dome  Thickness 
Sound  Velocity  in  the  fluid 
Density  of  the  fluid 
Operating  wavelength  in  the  fluid 
Plate  Modulus  of  dome  material 
Density  of  dome  material 


SYMBOL  dimension"^ 

a  L 

b  L 


h 

c 

0 


X  L 

+  -1  -2 

Eg  ML  -^T  ^ 

Os 


Dimensions  are  given  in  the  mass,  length,  time  system  with 
dimensions  M,  L,  T  respectively. 

4-  F 

E  =  - ;  E  is  Young’s  Modulus,  y  is  Possion's  Ratio 


In  order  to  cover  a  practical  range  of  sonar  applications, 
it  is  necessary  to  evaluate  the  various  performance  criteria  for 
a  wide  range  of  values  for  the  dome -transducer  parameters.  The 
direct  presentation  of  this  data  would  require  considerable  space 
and  be  quite  cumbersome  to  apply  directly  to  a  design  problem. 
Furthermore,  trends  in  the  data  would  be  difficult  to  recognize. 
The  problem  of  data  presentation  is  simplified  considerably  by 
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grouping  the  dome -transducer  parameters  in  dimensionless  products. 

Performance  criteria  then  can  be  displayed  as  functions  of  the 

dimensionless  products.  This  procedure  enables  the  designer  to 

recognize  trends  in  the  data,  and  the  volume  of  data  is  reduced 

to  a  set  of  compact  curves  which  can  be  used  efficiently.  Apply- 

2 

ing  the  methods  of  dimensional  analysis,  using  some  insight, 
and  performing  a  few  trial  and  error  calculations  leads  one  to 
adopt  the  following  dimensionless  products: 

1.  h/x,  dome  thickness  product 

2.  (b-a)/x,  dome -transducer  spacing  product 

2 

3.  E  / pc  ,  dome  stiffness  product 

s 

4.  p  /p,  dome  density  product. 

s 

The  product  2a/X,  the  ratio  of  transducer  diameter  to 
wavelength,  and  the  above  products  form  a  complete  set  of 
dimensionless  products.  Typical  values  of  2a/X  for  current 
sonar  systems  range  between  8  and  12.  A  nominal  value  of  11 
was  used  for  the  computations  included  in  this  technical  memo- 

2  a 

randum.  Results  for  other  values  of  t —  in  the  range  8-12  will 

A 

be  similar  qualitatively  to  the  results  presented  here. 


The  word  ’’product"  will  be  used  to  denote  dimensionless  groupings 
of  the  dimensioned  system  parameters  (h,  X,  etc.) . 
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III.  RESULTS 

This  section  contains  numerical  results  which  demonstrate 
the  dependence  of  performance  criteria  on  the  dimensionless 
products  of  the  dome -transducer  model.  Linear  plots,  most  of 
V7hich  have  the  thickness  product  h/x  as  the  independent  variable, 
are  employed  to  present  the  performance  criteria  data.  Each  plot 
contains  a  family  of  curves  generated  by  varying  one  of  the 
dimensionless  products  while  the  other  two  products  remain 
constant.  The  range  of  values  of  the  dimensionless  products  is 
chosen  to  span  the  areas  of  current  interest  and  demonstrate  the 
behavior  of  each  performance  criterion.  The  ranges  of  the 
dimensionless  products  are 


DIMENSIONLESS  PRODUCT 


RANGE  OF  VALUES 


h/x 

0.0 

(b-a)/x 

1.5, 

Pg/p 

2.0, 

E,/pc2 

90,  90  X 

,-2 


.-4 


The  thickness  product,  h/x,  spans  the  practical  limits  of  metal 
and,  possibly,  plastic  domes.  The  range  of  (b-a)/x  brackets  the 
minimum  and  maximum  dome-transducer  separation  of  the  AN/SQS-4, 
AN/SQS-23,  and  AN/SQS-26  systems.  Values  of  the  density  product, 

p  /o,  include  steel,  plastic,  and  some  wire-reinforced  rubber 

^  2 
compounds.  The  stiffness  product,  E  / oc  ,  has  a  range  which 

s 

covers  most  practical  metals  and  some  plastics.  Approximate 
values  of  the  dimensionless  products  for  a  system  such  as  the 
AN/SQS-26  are  h/x  =  0.02,  (b-a)/x  =  2.0,  p^/p  =  7.9,  and  E  /pc^  = 
90  (metal  dome) . 

Typically,  a  set  of  three  plots  is  shown  for  each  perform¬ 
ance  criterion.  In  the  first  plot,  density  product  p  / o  and 

s 
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stiffness  product  E  / pc“  are  constant;  spacing  product  (b-a)/x 

is  varied  to  produce  the  family  of  curves.  Values  of  density 

product  p  / p  are  shown  in  the  second  plot  for  constant  values 
^  2 

of  stiffness  product  E  /pc  and  spacing  product  (b-a)/X.  The 

s 

third  plot  is  generated  with  different  values  of  stiffness 
2 

product  E  /pc  for  constant  p  /p  and  (b-a)/X. 
s  s 

A.  Directivity  Index 

Directivity  Index  is  displayed  in  Figs.  2,  3,  and  4. 

Data  in  Fig.  2  indicate  decreasing  directivity  for  increased 

values  of  the  spacing  product  (b-a)/x  and  the  thickness  product 

h/x.  The  effect  of  the  (b-a)/X  is  less  pronounced  than  the 

effect  produced  by  h/x.  Figure  3  shows  the  improved  directivity 

produced  by  lower  values  of  the  density  product  o  /p.  As  shown 

s  2 

in  Fig.  4,  the  present  range  of  the  stiffness  product  E  /pc  has 
no  effect  on  Directivity  Index. 


B.  Major  Lobe  Width 

Results  for  Major  Lobe  Width  are  shown  in  Figs.  5,  6, 
and  7.  Major  Lobe  Width  increases  as  a  function  of  the  thickness 
product  h/x  but  decreases  as  a  function  of  the  spacing  product 
(b-a)/x  for  the  large  density  product  p  /p  as  shown  in  Fig.  5. 

o 

However,  for  small  p  /p  no  dependence  is  seen  on  either  h/x  or 

O 

(b-a)/x  as  shown  in  Fig.  6.  Results  in  Fig.  7  indicate  that 

2 

Major  Lobe  Width  is  independent  of  the  stiffness  product  E  / pc  . 

s 

C.  Side  Lobe  Level 


Data  for  Side  Lobe  Level  are  given  in  Figs.  8,  9,  and 

10.  The  effects  of  the  spacing  product  (b-a)/x  and  the  thickness 

product  h/x  for  the  constant  density  product  p  /p  and  the  stiff- 

2  ^ 

ness  product  E  /pc  are  shown  in  Fig.  8.  The  level  of  the  side 

s 

lobes  increases  for  increasing  values  of  h/x  and  (b-a)/x.  From 
Fig.  9,  reduced  values  of  density  product,  o  /p,  appear  to  produce 
side  lobes  which  are  lower  than  those  of  a  bare  transducer. 
However,  as  shown  in  the  farfield  directivity  patterns  in  Figs. 
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11  and  12,  this  phenomenon  is  associated  with  the  presence  of 
two,  nearly  equal,  side  lobes,  and  what  appears  as  an  improve¬ 
ment,  actually  compounds  the  problem.  Doubling  of  the  side  lobes 
appears  to  occur  for  values  of  (b-a)/x  of  1.5  -  2.5  .  Figure  10 

adds  data  indicating  the  independence  of  Side  Lobe  Level  on  the 

2 

stiffness  product  E  / oc  . 

s 

D .  Dome  Transmission  Coefficient 

Figures  13,  14,  and  15  display  results  for  Dome  Trans¬ 
mission  Coefficient.  Figures  13  and  15  show  transmission  to 

be  independent  of  the  spacing  product  (b-a)/x  and  the  stiffness 
2 

product  Eg/ pc  .  Strong  dependence  on  the  density  product  is 
sho^rm  in  Fig.  14.  For  values  of  the  density  product  i  /p 
approaching  unity,  transmission  becomes  independent  of  all 
products,  as  expected. 

E .  Dome  Transfer  Coefficient 

In  Figs.  16,  17,  and  18,  Dome  Transfer  Coefficient 
is  plotted  with  the  spacing  product  (b-a)/x  as  the  independent 
variable.  The  data  show  Dome  Transfer  Coefficient  to  be  a  damped, 
periodic  function  of  (b-a)/x,  with  period  x/2.  This  behavior 
is  attributed  to  strong  standing  waves  set  up  in  the  dome- 
transducer  model.  The  standing  wave  patterns  in  real  domes 
should  not  be  so  pronounced;  however,  such  geometry  effects 
cannot  be  completely  ignored.  Analytical  models  with  more 
realistic  geometry  will  yield  better  determiLnations  of  the 
effects  of  system  geometry  on  performance. 

F .  Nearfield  Roughness 

Results  for  Nearfield  Roughness  are  presented  in 

Figs.  19,  20,  and  21.  Figures  19  and  21  show  Nearfield  Roughness 

to  be  independent  of  the  spacing  product  (b-a)/x  and  the  stiffness 
2 

product  E  / pc  .  In  Fig.  20  an  increase  in  Nearfield  Roughness 
is  noted  for  Pg/p  >  2.0  .  For  smaller  values  of  the  density 

product  p  / p  Nearfield  Roughness  appears  to  be  independent  of 

o 

all  dimensionless  products. 
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FIXED  DENSITY  PRODUCT  ^ 
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IV.  CONCLUSIONS 

The  numerical  results  presented  in  this  memorandum  provide 
a  quantitative  measure  of  the  dependence  of  certain  criteria  of 
transducer  performance  on  the  geometrical  and  material  parameters 
of  the  dome -transducer  model.  While  these  results  do  not  apply 
directly  to  existing  sonar  equipments  because  of  the  unrealistic 
dome  geometry,  some  qualitative  conclusions  can  be  made  in  regard 
to  the  effects  of  the  sonar  dome  on  system  performance. 

In  general,  the  dome  causes  a  degradation  of  all  the  assumed 

dome  performance  criteria,  as  compared  to  a  transducer  without 

dome.  The  results  indicate  that  the  strongest  dependence  of  the 

various  performance  criteria  is  on  the  dome  thickness  product 

h/x  amd  dome  density  product  o^/o.  Weaker  dependence  on  the 

dome  spacing  product  (b-a)/X  is  seen,  while  the  stiffness  pro- 
2 

duct  Eg/ pc  exerts  little  influence  on  the  performance  criteria. 

It  is  helpful  in  explaining  the  results  to  consider  the 
acoustic  impedance  of  the  dome.  The  point  acoustic  impedance  is 
defined  as  the  ratio  of  net  acoustic  pressure  to  normal  particle 
velocity  at  the  mean  surface  of  the  dome.  The  point  impedance 
of  each  mode  of  the  dome  is  given  by 


^n  = 


■iiDo  h  + 
s 


ih 

11)  ^ 


iE  h-^  / 

S  r_4 

- T  [n 

12u)b 


2n^  +  1]  + 


ihn^Eg^ 

u)b^(p  tu^b^  -  n^E  ) 
'  s  s 


This  expression  is  obtained  by  solving  for  the  ratio  of  net  acoustic 
pressure  to  normal  particle  velocity  in  the  shell  equations  des¬ 
cribing  the  motion  of  the  dome.  With  some  manipulation  this 
expression  is  similar  to  the  result  obtained  by  Heckl^  and  later 
employed  by  Barger^  with  the  exception  that  Flugge's  shell  equations 
rather  than  Kennard’s  have  been  used  in  this  study.  (See  Ref.  5.) 
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From  this  expression  it  is  obvious  that  the  point  impedance  is 
directly  proportional  to  dome  density  and  thickness  h.  Point 

impedance  is  inversely  proportional  to  the  dome  radius  b;  however,  \ 

the  dependence  of  on  b  is  not  a  strong  one.  It  is  seen  in  the 

results  that  each  of  the  performance  criteria  degrade  as  Z^  is 

increased  (increasing  h/>.  or  0^/0).  An  increase  in  (b-a)/\ 

usually  leads  to  a  slight  improvement  in  the  performance  criteria. 

Based  on  these  results  it  appears  that  the  most  important 

dimensionless  parameters  in  the  dome -transducer  system  are  the 

thickness  product  h/x  and  the  density  product  p  /p,  and  that 

the  initial  guidelines  could  be  compacted  further  by  computing 

the  dependence  of  the  performance  criteria  on  a  new  dimensionless 

variable  obtained  as  the  product  h/x  •  0^/0  .  From  the  relation 

Xi  =  2nc,  where  lu  is  the  angular  frequency  and  c  is  the  speed 

of  sound  in  water,  the  new  dimensionless  variable  6  can  be 
(UP  h 

obtained  as  6  “  .  This  expression  is  seen  to  be  approxi¬ 

mately  Z  normalized  to  the  specific  impedance  of  water,  or  the 
impedance  "mismatch"  of  the  dome . 

The  results  contained  in  this  memorandum  support  the 
following  specific  conclusions  regarding  the  effects  of  the  dome 
on  transducer  performance. 

1.  The  performance  criteria  of  the  dome- transducer  are 

a  strong  function  of  the  variable  '.ijD^h/pc,  an  approximate  measure 
of  thc‘  point  impedance  of  the  dome.  Improved  acoustic  performance 
is  obtained  by  using  dome  materials  having  a  surface  density 
p  h  as  small  as  possible.  This  conclusion  is  intuitively  obvious 
for  such  criteria  as  Dome  Transmission  Coefficient,  but,  more 
important,  is  valid  for  farfield  distribution  characteristics 
typified  by  Directivity  Index. 

2.  Performance  criteria  are  not  a  strong  function  of 
dome- transducer  spacing  but  improve  slightly  with  increased 
spacing.  Variation  in  power  transferred  to  the  water  as  a  result 
of  reflections  from  the  dome  are  undoubtedly  exaggerated  by  the 
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dome  geometry  used  in  this  study.  However,  some  interference 
can  be  expected,  particularly  in  domes  having  broad,  near- 
concentric  windows  such  as  the  AN/SQS-26  dome.  It  should  be 
noted  that  scattering  from  a  periodic  structure  such  as  a  typical 
dome  supporting  structure,  not  considered  in  this  memorandum, 
can  complicate  further  the  effects  of  reflections. 

3.  Performance  criteria  are  not  dependent  on  the  elastic 
properties  of  the  dome  material.  Hence,  high-strength  materials 
should  not  degrade  acoustic  performance.  It  can  be  concluded 
that  a  low-density,  thin  dome  made  of  high-strength  material 
should  result  in  improved  acoustic  performance  and  satisfy 
strength  and  durability  requirem.ents . 

The  sets  of  curves  contained  in  this  memorandum  are  repre¬ 
sentative  of  the  information  necessary  for  the  design  of  sonar 
domes  on  a  scientific  basis.  Improved  design  guidelines  await 
the  development  of  more  sophisticated  analytical  models. 
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APPEtJDIX  A.  THE  SOUND  PRESSURE  FIELD  PRODUCED  BY  A  TRANSDUCER 

RADIATING  IN  A  CONCENTRIC  DOME 


The  mathematical  model  illustrated  in  Figure  1  consists 
of  an  inf ini te "length  cylinder  and  concentric  shell.  A  velocity 
distribution  independent  of  the  axial  direction  is  specified  on 
the  surface  of  the  cylinder,  and  it  is  desired  to  find  the 
resulting  sound  pressure  field  in  the  fluid  in  Regions  1  and  2 
of  Figure  1.  The  boundary  value  problem  as  stated  is  a  three 
region  problem  (two  fluid  regions  and  the  shell).  In  general 
partial  differential  equations  must  be  solved  both  in  the  two 
fluid  regions  and  in  the  region  of  the  shell  material,  and 
boundary  conditions  must  be  matched  at  shell-fluid  interfaces 
as  well  as  at  the  cylinder-fluid  interface.  However,  the 
problem  is  considerably  simplified  by  using  the  results  of 
elastic  shell  theory  to  describe  the  motion  of  the  shell's  middle 
surface,  thus  eliminating  the  necessity  of  solving  the  full 
partial  differential  equations  of  linear  elasticity  for  the  shell 
material.  Further,  since  the  shell  is  thin  relative  to  other 
gecxnetrical  parameters, the  boundary  conditions  at  shell-fluid 
interfaces  can  be  matched  at  the  middle  surface  of  the  shell. 

A.  Sound  Pressure  Field  in  the  Fluid 

It  is  assumed  that  in  Regions  1  and  2  the  acoustic 
pressure  p  obeys  the  scalar  wave  equation 


1  > 
r  ar 


(r 


+  JL  ife.  _  L  i!2 

*  ^  r\  r\ 


(1) 


where  r  and  0  are  the  space  coordinates! t  is  the  time,  and  c  is 
the  propagation  velocity  in  the  fluid. 

For  a  harmonic  time  dependence  the  general  solution 
to  Eq.  1  is 
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FIG.  I  “  GEOMETRY  FOR  A  CYLINDRICAL  TRANSDUCER 
RADIATING  IN  A  CONCENTRIC  DOME 
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to 

V  =  y  [A„  H^^^(kr)  +  cos  n0  (?) 

il  Xl  li  Xl 

n=0 

where  and  are  the  nth  order  Hankel  functions  of  the 

n  n 

first  and  second  kind,  respectively,  u)  is  the  angular  frequency, 
k  =  uj/c,  and  A  and  B  are  constants  which  will  be  determined 
by  boundary  conditions.  Equation  2  involves  only  cosine  terms 
as  it  is  anticipated  that  problems  to  be  analyzed  will  possess 
symmetry  about  0=0.  In  Region  1  the  sound  pressure  field  is* 


Pi  “  II  H^®^(kr)]  cos  ne 

n=0 


(3) 


The  sound  pressure  and  the  radial  component  v^,  of  the  fluid  velocity, 
with  harmonic  time  dependence,  are  related  by 


_i!r 

p  3r  dt 


=  -iujv^ 


(^) 


At  the  shell -fluid  interfaces  this  component  of  the  fluid  velocity 
must  equal  the  velocity  of  the  shell.  (This  restriction  does 
not  apply  to  velocity  components  tangential  to  the  interface 
because  the  fluid  is  regarded  as  inviscid  and  hence  will  not 
support  shear  forces.)  Denoting  the  velocity  on  the  cylinder 
by  V^.  and  the  radial  component  of  the  velocity  of  the  shell  by  w, 
the  interface  conditions  at  r=a  and  r=b  are 
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(5) 


(6) 


*The  time  component  will  be  suppressed  in  the  following. 
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In  the  case  of  Interest,  a  section  of  the  cylinder  of 
angular  extent  2cp  is  phased  to  represent  a  plane  velocity  dis¬ 
tribution  of  uniform  amplitude,  while  the  remainder  of  the 
cylinder  is  rigid.  This  condition  can  be  written  as 

v,  =  ^=os  e  »  ,  iei<,  . 

-  0  .  lel>tp  .  (7) 

4  O  ^  A  A 

A  Fourier  series  for  the  function  cos  0  e  may  be  obtained 

by  differentiating  the  identity, 

CD 

^Ikrcos  e  ,  y  .  1"  J  (kr)  cos  ne 

A — I  n  n 

n=0 

where  e^=l,  c=2  for  n  >  o,  and  J  is  the  Bessel  function.  Then 
o  n  '  n 

OB 

^r  “  ^  ®n  Jn'(ka)  cos  n0  ,  10l<cp  , 

n=0 


=  0 


10l>cp  ,  (8) 


where  the  prime  indicates  a  derivative  with  respect  to  the 
argument .  Equation  8  can  be  represented  by  the  Fourier  series 


where 


and 


ipc  Z_j 
n=0 


6^  cos  n0 


y  I*"  J  '  (ka) 
n  n  /  ■  m  m  '  ' 


m=0 


(9) 
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jnj  ,  n/o  ;  M 


The  substitution  of  Equations  2  and  9  into  Eq.  5  results  in 

00  eo 

y  [A  H^^^'(ka)  +  B  '  (ka)]cos  nG  =  /  6_  cos  nG .  (10) 


k 

imp 


n=0 


n=0 


Noting  that  the  cosine  functions  are  linearly  independent  and 
that  k/u)  =  l/c,  Eq .  10  implies 

An  H^^^'(ka)  +  ^  ^ 

The  radial  component  of  the  shell  velocity  can  be  represented  in 
the  form 

00 

w  =  ^  cos  nG  .  (12) 

n=0 

The  boundary  condition  given  by  Eq.  6  then  becomes 
00  00 

y-  -  y  [A„  H'5^^'(kb)  +  B„  '  (kb)]cos  n0  =  y  G„  cos  n0  (13) 

iiup  ^  ‘■nn  nn  ^  j  n  v 

n=0  n=0 

or 

^n  ^n^^ '  K  ”n^^  '  ^n 

since  the  cosine  functions  are  linearly  independent. 
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In  writing  the  general  solution  for  the  sound  pressure 
Pg  outside  the  shell  it  is  necessary  to  account  for  the  fact  that, 
for  large  r,  the  pressure  field  should  have  the  form  of  an  out¬ 
going  wave.  Since 


Llm  H(l)(kr)|  =  ishVt 

represents  an  outgoing  wave,  the  proper  form  for  Pg  is 

00 

Pg  “  ^  H^^^(kr)  cos  n6  .  (15) 

n=0 

The  interface  condition  on  the  velocities  of  the  fluid  and  shell 
on  the  outside  surface  of  the  shell  will  be  matched  at  the  middle 
surface  of  the  shell,  as  in  Eq.  6.  Hence 


1  “Pa  1 

iujp  dr  ' 

00 

00 

—  7 

iiup  ^ 

n0  =  7  f'n  ’ 

(16) 

n=0 

3 

II 

O 

and,  noting  the  linear  independence  of  the  cosine  functions, 

Cn  '  (kb)  =  icp  .  (17) 

Equations  11,  14,  and  17  give  three  sets  conditions  on 
the  four  sets  of  unknowns,  A^,  B^,  and  A  consideration 

of  the  motion  of  the  shell  is  necessary  to  obtain  the  requisite 
fourth  condition. 
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B .  Motion  of  the  Shell 


Let  the  radial  and  tangential  displacements  from 
equilibrium  of  the  middle  surface  of  the  shell  be  denoted  by 
Uj,  and  Ug,  respectively.  Flligge's  shell  equations  are 


^2 


*^9  .  »“r, 


(18) 


at2 


p-  '  p+  .  ^s  A 

h  ^  *-30 


p-  k2  ^4  ^2 

8_h  a  u  3  u„ 

+  u  ]  -  -2-if  [ - 1  +  2  — ^  +  u  1  (19) 

12b^  30^  30 


where  h  is  the  thickness  of  the  shell  and  b  is  the  radius  of  its 
middle  surface,  p  is  the  mass  density  of  the  shell  material  and 
E  =  E/ (1-v  ),  E  being  Young's  modulus  and  v  Poisson's  ratio. 

O 

p_  and  p^  are  the  pressures  on  the  inside  and  outside  shell 
surfaces,  respectively,  so  p_  =  and  p^  =  [Pglr^b* 

Since  w  is  the  radial  velocity  of  the  shell, 

3u 

w  =  =  -imUj.  .  (20) 


The  tangential  velocity  v  of  a  shell  element  is  represented  as 

3Up 

=  aF"  =  ■^'"“0  (21) 

For  symmetry  about  0=0,  the  Fourier  series  representation  for  v 
is  of  the  form 


V 


I 


Dn  Sin  n0  . 


The  shell  Equations  18  and  19  can  be  differentiated  with 
respect  to  time  to  give  conditions  on  u  and  v. 


(22) 
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-U)  p  V  =  -p 


2 

-'D  p 


b2  Lae 


+  w] 


(23) 


ts  A  ^  2  ^  +  w]  .  (24) 

i2b^  ze  hQ 


If  Equations  4,  16  and  the  Fourier  series  for  w  and  v  are 
substituted  into  Equations  23  and  24  and  the  linear  independence 
of  the  trigonometric  functions  employed,  the  following  conditions 
are  obtained. 


-iiu  [Hj![^^(kb)  +  H^^\kb)Bj^]+im  H^^^(kb)  C 


h  E_  n 


- ^ - 

n  ^2  n 


+  h  p,  -  %  [1  +-^(n'^-2n^+l)] 

I  ®  b  ia> 


Ps  ■  ^  "An  -  ^  =  0 


Gn  =  0  (25) 


(26) 


C .  Determination  of  the  Coefficients 


Equations  13,  15,  18,  25,  and  26  represent  a  set  of 

simultaneous  algebraic  equations  in  the  unknowns  A  ,  B  ,  C  ,  D  , 

n  11  n  n' 

and  This  set  of  equations  is  readily  solved  to  complete  the 

solutions  for  the  sound  pressure  field  in  Regions  1  and  2  given 
by  Equations  3  and  15  respectively. 


As  an  example  of  the  form  of  the  coefficients,  the  solution 
for  the  set  of  coefficients  which  are  associated  with  Pg 
(Eq.  15)  is 
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C_  = 


n 


n 


Hj;'(ka)  +  h  7^  H^'(kb)  [J^^Ckb)  N^(ka)  -  Nj;(kb)  J^(ka)] 


(27) 


where 


nb 


^  2 
2c  p 


2  4  2 

U)  Ps'7^tl  +  7~f  2  -  2n  +  1)] 


12b  ‘ 


b^(b^u)%g  -  Eg  n^) 


This  form  is  typical  of  the  remaining  coefficients.  It  is  of 

interest  to  note  that  for  h  =  0,  the  C  reduce  to  the  coefficients 

n 

obtained  for  a  transducer  radiating  without  a  dome. 
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APPENDIX  B.  DETERMINATION  OF  TRANSMISSION  COEFFICIENT 

It  is  often  helpful  to  relate  certain  phenomena  associated 
with  sound  transmission  through  shells  to  their  analogues  in  the 
theory  of  plane  wave  transmission  through  plates.  In  the  latter 
area  of  study  there  are  basically  two  t3rpes  of  problems,  one 
in  which  the  velocity  distribution  is  specified  on  the  surface 
radiating  the  pressure  wave  and  the  other  in  which  the  form  of 
the  plane  wave  incident  on  the  plate  is  specified.  Of  course, 
the  first  type  of  problem  more  nearly  approximates  reality  in 
that  it  takes  into  account  coupling  between  the  plate  and  the 
radiating  surface  (multiple  reflections,  etc.).  But  the  second 
type  of  problem  is  particularly  well  suited  for  studying  proper¬ 
ties  of  the  plate  itself  and  is  usually  the  problem  seen  in  the 
literature . 

In  the  second  type  of  problem,  the  form  of  the  incident 
wave  is  given  and  the  waves  transmitted  through  and  reflected 
from  the  plate  are  calculated.  Then  a  plate  transmission  co¬ 
efficient  T  can  be  defined  as  the  ratio  of  the  amplitude  of  the 
transmitted  wave  to  the  amplitude  of  the  incident  wave.  If  the 
incident  wave  is  not  a  plane  wave  but  is  a  type  of  wave  which 
would  be  emitted  by  a  finite  source,  the  transmission  coefficient 
can  be  defined  as  the  ratio  of  power  transmitted  through  the 
plate  to  power  incident  on  the  plate.  In  this  case  the  trans¬ 
mission  coefficivent  is  a  function  of  the  form,  but  not  of  the 
amplitude,  of  the  incident  wave.  For  a  given  form  of  incident 
wave  the  transmission  coefficient  depends  only  on  plate  param.eter 

The  basic  problem  described  in  Ref.  1  is  analogous  to  the 
first  type  of  plane  wave  problem  mentioned  above  in  that  the 
coupling  between  the  cylindrical  shell  and  the  cylindrical  radi¬ 
ating  surface  transducer  is  included.  However;  in  this  problem 
it  is  impossible  to  separate  the  resultant  field  between  the 
transducer  and  shell  into  incident  and  reflected  parts  and  thus 
define  a  transmission  coefficient.  Nevertheless,  it  is  possible 
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to  formulate  a  problem  analogous  to  the  second  type  of  plane  wave 

problem  mentioned  above  by  specifying  the  form  of  the  incident 

cylindrical  wave  and  calculating  the  transmitted  wave.  Then  a 

transmission  coefficient  can  be  defined  as  the  ratio  of  power  ' 

transmitted  through  the  shell  to  power  incident  on  the  shell. 

Of  course,  in  this  type  of  problem  the  wave  reflected  inward  from 
the  shell  is  not  scattered  back  toward  the  shell;  that  is, 
multiple  reflections  are  not  included.  The  transmission  coeffici¬ 
ent  will  be  a  function  of  the  form  of  the  incident  wave*  however, 
it  should  not  vary  much  among  those  types  of  waveforms  which 
represent  directed  beams. 

A  typical  directed  beam  is  that  produced  by  a  cylindrical 
transducer  driven  in  the  manner  described  in  detail  in  Ref.  1. 

To  obtain  this  beam  the  elements  on  a  120°  section  of  the  cylin¬ 
der  surface  are  driven  with  the  phase  and  amplitude  distribution 

fcll 

of  a  plane  wave.  The  amplitude  and  phase  of  the  j  element  will 

be  denoted  by  S.  and  co. ,  respectively,  and  the  driving  and  response 

Ji  “idit 

will  be  assumed  to  have  an  e  time  dependence.  Tlien  the 

velocity  of  the  fluid  at  the  surface  of  the  transducer  can  be 
expanded  in  a  Fourier  series  in  the  0  coordinate  of  the  cylindri¬ 
cal  coordinate  system  r,  0,  the  surface  of  the  transducer  being 
r  =  a.  Thus 

03 

-iu)t  ^ 

■y  =  - -  6  cosn0  ,  (1) 

where  d  is  the  fluid  density  and 

CO 

6  =  imp—  sinn[ (2j -l)a+3 ]  -  sinn[ (2j -l)a-0]|  ,  (2) 

n  nn  .  j  L  J 

j=o 
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with  23  the  angular  width  of  the  elements  and  2p  their  angular 
spacing.  Tne  resulting  directed  beam  is  then 

09 

Pj,nc  "  Yj  \  cosnQ  ,  (3) 

n=o 

where  (hr)  is  a  Hankel  function  of  the  first  kind,  k  is  the 

ratio  of  lu  to  the  speed  of  sound  c  in  the  fluid,  and 

A„  =  '  (ka)  ,  (4) 

the  prime  denotj.ng  differentiation  with  respect  to  argument. 

In  order  to  calculate  a  transmission  coefficient  for  the 
shell  it  will  be  assumed  that  the  pressure  wave,  Eq.  (3),  can  be 
generated  within  the  shell  and  directed  at  the  shell  surface  so 
that  the  radiated  energy  will  be  partly  transmitted  through  the 
shell  and  partly  reflected.  The  transmitted  and  reflected  pressure 
fields,  p^  and  p^,  respectively,  have  the  forms 

09 

cosn0  ,  (5) 

n=o 


-i  M,t-  /■9^ 

p  =  e  y  B  H  '' (kr)  cosn9 

n=o 


(6) 


'2) 

where  H  •  ' (kr)  is  a  Hankel  function  of  the  second  kind,  and  B 
and  C  are  constants  which  will  be  determined  by  conditions  at 
the  shell,  at  r  =  b.  It  should  be  noted  that,  for  an  e~  Lime 

dependence,  the  Hankel  functions  represent  outgoing  waves 

and  the  H^^^  Hankel  functions  represent  ingoing  waves.  The 
pressure  field  inside  the  shell  is  the  sum  p.  +  p.- 
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The  method  of  analysis  proceeds  exactly  as  in  Ref.  1. 

The  pressure  fields  are  substituted  into  shell  equations  involv¬ 
ing  the  shell  velocity  components,  Fourier  series  are  ass'  .led  for 
these  velocity  components,  then  the  coefficients  ’.n  these  series 
and  in  the  series  of  Eqs .  (5)  and  (6)  can  be  determined. 


1 


D  A 
n  n 


L  D_ 


n 


(7) 


C 

n 


(8) 


In  these  expressions 


Dn  = 

4c  p 


E  ■•2  n  n 

-y  [IH - ^(n  -1)  ] — ^ 

12b^ 


2 

^s  ^ 


b^[b^m  0,-E  n^] 
^  s  s 


2 

h  is  the  shell  thickness,  p  its  mass  density  and  E  =  E/ (1-v  ), 
where  E  is  Young's  modulus  and  v  is  Poisson's  ratio.  The  pressure 
fields  are  completely  determined. 

The  average  power  inv ident  on  the  shell  and  transmitted 
through  the  shell  are  then 


Power  Incident  = 


mp 


''O'i 


Power  Transmitted 


=  L-  Y 

UJP  Z, 


n=o 


n 


> 
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where  s  ■ 
o 

is  defined 


1>  ^  n  >  1.  The  transmission  coefficient  T 

as 


n,  _  Power  Trfnsmitted 
Power  Incident 


(11) 
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